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In a double-well potential, a Bose-Einstein condensate exhibits Josephson oscillations or self- 
trapping, depending on its initial preparation and on the ratio of inter-particle interaction to inter- 
well tunneling. Here, we elucidate the role of the exchange symmetry for the dynamics with a 
mixture of two distinguishable species with identical physical properties, i.e. which are governed 
by an isospecific interaction and external potential. In the mean- field limit, the spatial population 
imbalance of the mixture can be described by the dynamics of a single species in an effective 
potential with modified properties or, equivalently, with an effective total particle number. The 
oscillation behavior can be tuned by populating the second species while maintaining the spatial 
population imbalance and all other parameters constant. In the corresponding many-body approach, 
the single-species description approximates the full counting statistics well also outside the realm 
of spin-coherent states. The method is extended to general Bose-Hubbard systems and to their 
classical mean- field limits, which suggests an effective single-species description of multicomponent 
Bose gases with weakly an- isospecific interactions. 



I. INTRODUCTION 

Multi-component Bose gases pQ exhibit a panoply of 
different quantum phases [2 -4 , which reflect the hierar- 
chy of inter- to intra-species interaction parameters that 
compete with a possibly species-dependent [5 external 
potential. The an-iso specificity, i.e. the distinct physical 
properties of the different species, is the main cause for 
the great increase of complexity that mixtures feature 
with respect to single-species Bose gases. 

A minimalist exemplary system that contains many of 
the building blocks for quantum many-body dynamics is 
the bosonic Josephson junction [6j[7], i.e. a Bose-Einstein 
condensate (BEC) in a double-well potential. Despite 
its simplicity, it features qualitatively different dynami- 
cal regimes - self-trapping and Josephson oscillations - 
already in the single-component case [6j I5HTT] , which has 
been realized experimentally with coupled spatial modes 
[12] [13] and hyperfine states [14] . When a bi-component 
BEC is loaded into the double- well [T5H2Q] , new features 
emerge: For example, chaos p~5j [21], phase separation 
transitions [22] as well as ferromagnetic behavior [23] can 
arise, relying on the distinct inter- and intra-species in- 
teractions in the system. 

If, on the other hand, all intra- and inter-species 
interactions are equal and all external potentials are 
species-independent, adding a second species to an ini- 
tially homo-specific system (or, equivalently, populating 
a second internal state) does, at first sight, not appear 
to change the system dynamics: The new species does 
not feature any distinct property in comparison to the 
present one. However, the mere assignment of a particle 
label - although invisible for the Hamiltonian - breaks 
the exchange symmetry of the quantum state of the sys- 
tem, with important consequences. 

Photons, for instance, do not interact, but indistin- 
guishable photons can interfere collectively. This leads 
to dramatic differences between the counting statistics 



of photons that can be distinguished by their polariza- 
tion or frequency (and thus do not interfere collectively) 
and identical photons [24H27] . The latter are governed 
by many-particle interference, such that, e.g. events with 
many particles in one mode are privileged [28] [29] . 

For the bosonic Josephson junction, the question arises 
how the combination of inter-particle interaction and 
a broken (or immaculate) exchange symmetry affects 
the dynamics. This investigation is the purpose of the 
present article: We compare a bi-species condensate with 
two species of identical properties, for which we pro- 
vide an experimental protocol, to a single-species con- 
densate. This reveals the consequences of broken ex- 
change symmetry in a paradigmatic interacting many- 
particle system. We show that the population of a sec- 
ond species can significantly affect the system dynam- 
ics: The tunneling behavior (self-trapping or oscillatory) 
can be switched by the manipulation of the species pop- 
ulations. This paves the road to the control of many- 
particle tunneling dynamics by merely populating differ- 
ent internal states. In the mean- field limit, the dynamics 
can be understood from the broken phase-coherence be- 
tween the wells, and an exact single-species description 
is formulated. In the fully second-quantized formulation, 
the bi-species Fock-state counting statistics is efficiently 
approximated by the single-species model. We also in- 
vestigate the relaxation of isospecificity, i.e. species with 
slightly different parameters, and show that the single- 
species description still captures the essential dynamics 
within a wide range of parameters. The generalization to 
larger systems therefore suggests an efficient description 
of weakly an-isospecific multicomponent Bose gases by a 
single species. 

We first review the main properties of the single-species 
bosonic Josephson junction - in the discrete two-mode 
approach - in Section [TTJ The physical consequences of 
the addition of a second species can be absorbed by an 
effective particle number or by an effective tunneling cou- 
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FIG. 1: (color online) Phase-space of the bosonic Josephson junction in the z^ — ^-plane, for different values of N^U/Q. 
For convenience, we show the range —tt/2 < 9^ < 3/2tt. Thick solid lines denote the separatrices between the regions of 
qualitatively distinct behavior, dotted lines show exemplary trajectories, and fix points are denoted by red crosses. The fix 
point at (z^, 9^) = (0, tt) (N^U/Q < 1) becomes unstable for N^U/Q > 1. In this regime, two localized fixed points with z^ ^ 
emerge. Solutions in the blue areas are characterized by (z^) = (0^) = 0, solutions in the white areas lead to (z^) = 0, (0^) = n. 
Self-trapped solutions are either 7r-modes (orange area, (z^) ^ 0, (9^) = 7r), or phase-running modes (green area, (z^) ^ 0, 
unbounded phase 9^). 



pling, as shown in Section III A complementary quantum 
calculation in the two-mode Bose-Hubbard model is per- 
formed for spin-coherent and for Fock-states in Section 
|IV[ which confirms the validity of the approach beyond 
the mean- field limit. The generalization of the results 
to a general Bose-Hubbard Hamiltonian and to its mean 
field limit is presented in Section [V] before we conclude 
in Section [VH 



is the dimensionless quotient N^U/Q, and we will mea- 
sure all frequencies in units of Q. Since the scaling of Q 
will also be discussed, we keep Q as an explicit parame- 
ter. Eq. represents a self-trapping equation [3Tj [32], 
which has been studied extensively [6] [TQl [TTJ [15] [30] f33l 
135] . We review here the essential properties that will be 
useful for our analysis of the two-component case; in par- 
ticular, we discuss a non-linear scaling property for the 
population imbalance between the wells. 



II. SINGLE-SPECIES GROSS-PITAEVSKII 
EQUATION AND NON-LINEAR SCALING 

In order to obtain analytic insight into the problem, 
we follow a two-mode treatment of the double-well po- 
tential, which is well-justified for very low temperatures 
and small condensate depletion [30 , and which was quan- 
titatively verified experimentally in Ref. [13]. The dy- 
namics of the macroscopic wavefunction that describes 
a (single-species) condensate in the symmetric double- 
well potential is then governed by the discrete Gross- 
Pitaevskii equation [9], [10] (we set h = 1), 



(1) 



where j = 1,2; ipi (i/j 2 ) is the amplitude of the wave- 
function in the left (right) well, such that \ipj\ 2 is the 
expectation value of the number of particles in the jth 
well, and = \ipi\ 2 + I V^l 2 is the conserved total num- 
ber of particles. The tunneling coupling Q and the in- 
teraction U can be inferred from the overlap integrals 
of the exact localized wavefunctions of the double-well 
potential [10 . The relevant parameter for the dynamics 



A. Bloch-sphere description 

The dynamics of the system can be re-formulated as 
the motion of a vector on a Bloch-sphere of radius 
[T5] . where 
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i.e. the z^-component of is related to the population 
imbalance between the wells, while the relative phase be- 
tween the left and right component of the condensate, 



64 = arg(^ 2 ), 



(3) 



is encoded in 



z^sinfy. (4) 
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The equation of motion Eq. ([I]) assumes the form of an 
optical Bloch equation with non-linear terms [20] , 



dt 



z^ + U z^ x^ 
-ft 



and the conserved energy reads 
U 



(5) 



(6) 



Depending on the ratio N^U/Vt, the phase space of 
the system assumes a distinct structure [10]. The three 
topologically distinct cases are shown in Fig. [T] For 
N^U/Vt < 1 (left panel), in the Rabi regime, all solu- 
tions are oscillatory with (z^) t = 0. There are two el- 
liptic fix points with z^ = and x^ = ±A/^, i.e. the 
average value of the phase 6^ is either or n. In 
the parameter range 1 < N^U/Vt < 2 (middle panel), 
the fix point (x^,y^,z^) = (A^,0,0) remains stable, 
while (x^jy^jZ^) = (— iV^,0, 0) becomes unstable. Two 
new fix points emerge, which are related to trapped n- 
mode oscillations with (z^) t ^ and (0^)t = tt- For 
N^pU/Q > 2 (right panel), trapped phase-running modes 
emerge, for which the phase is unbounded. When the in- 
teraction is increased further, these phase-running modes 
dominate the picture more and more, whereas the oscil- 
lation region shrinks - the Fock regime is attained. 



B. Second-order differential equation 

The overall behavior of the system is governed by the 
phase-space structures shown in Fig. [I] In general, the 
self-trapping or oscillatory behavior of a solution depends 
on the initial preparation of the system and on the ratio 
N^U/Vt. If one focuses on the population imbalance z^ 
only, a scaling symmetry relates solutions with different 
initial conditions and different values of U and Q to each 
other. For this purpose, we re- formulate the dynamics of 
z^ by taking the time-derivative of Eq. ^ [15] , 



-z^(ft 2 -2 UH)-U 



z 3 



(7) 



where H is the value of the conserved energy, Eq. mj. 
In terms of the initial preparation of z^ and i^, Eq. (m) 
reads 



(8) 



where x^(0) can be expressed through 2^(0) and i^(0), 
except for its sign: 



a V ,(0) = ±\/^-^(0)2-^ 



(9) 



Since the transformation U —> —U is equivalent to 
x^ —x^ [10], we choose here the convention x^(0) > 
(positive square root) without loss of generality. This im- 
plies that the initial phase 0^ lies between — tt/2 and 7r/2, 
see Eq. Q. Any self-trapping then comes along with a 
phase-running mode. 



C. Non-linear scaling of population imbalance 

From the second-order differential equation of motion 
Eq. ([sj) a scaling symmetry for the solutions can be de- 
rived by merely redefining the unit of time, 



2^(iV<0, 2^(0), ity(0), U, ft, t) 

,^ ^(°) u ft ^ 

= ^(^,^(0), -^7—, — , : ,at), 



a a a 



(10) 



for any positive scaling parameter a. Also the particle 
number can be scaled, 

z,f, (N^ , (0) , z^ (0) , U, ft, t) 

= -2:^(7^, 7^(0), 7i^(0), —,^,t), (11) 
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when the interaction strength U is scaled reciprocally 
with iV-0. Solutions on Bloch-spheres of different radii 
are thus related to each other. 

For vanishing interaction, [7 = 0, Eq. ([8| becomes a 
linear differential equation, and 



MA^,Mo),o,*y = o,ft,t) 

= ^(A^, 0*0(0), 0, U = 0, SI, t), 



(12) 



where < /3 < N^/z^(0) is a real parameter that scales 
the initial population imbalance, keeping the Bloch- 
sphere radius constant. Rabi oscillations of different am- 
plitude are related to each other, and the scaling Eq. ( 12 ) 



reflects the trivial topology of the phase space for U = 0, 
where all solutions z^(t) are sinusoidal. 

In general, a finite non-linearity [/ / breaks the 
linear scaling property Eq. (12): Different initial con- 



ditions for the population imbalance 2^(0) lead to quali- 
tatively different solutions, as already apparent from the 
phase-space structure for N^U/Vt > 1 in Fig. [I] Eq. (J8|, 
however, allows to relate solutions with different initial 
population imbalance to each other. More precisely, for 
vanishing initial derivative, i(0) = 0, a solution with 
2^,(0) = zq can be related to a scaled solution with 
2-0(0) = Co, such that 

2:0(^0, zo,0,U,Sl,t) = ^0(AT0,Co,O,C/,O,£), (13) 
Co 

when the interaction U and the tunneling coupling O are 
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scaled according to 



U 



= u 



(14a) 



u 2 4 

4 Co 2 



a) 



u 

' ~2 



Co 2 , 



(14b) 



remember that we assume 0^(0) = 0, or, equivalently, 
x^(0) > 0,^(0) = 0. A change of the initial population 
imbalance z^(0) from zo to Co can thus be compensated 
by choosing a scaled interaction U and tunneling CL. The 
relative phases of the two solutions, 0^, or, equivalently, 
the motion of and are, however, not linearly related 
to each other. Therefore, this scaling property is not 
apparent from the phase-space picture, since it relates 
solutions with possibly different topology to each other, 
and U/Cl ^ U/Cl, in general. Note that Cl may become 
complex. 



The scaling of the interaction [/, explicit in Eq. ( 14a 
can be understood rather intuitively: When we relate 
a solution with, say, z^(0) = zo = (all particles in 
the left well) to a solution with z^(0) = Co < A/^, the 
effect of the interaction U will be weaker for the system 
with smaller initial imbalance Co • To retrieve the original 
behavior of the solution with z(0) = A^, this effect has 
to be counter-balanced by the upscaling of U. 

The scaling of the frequency Cl is more intricate. In 
the first place, Cl ^ Cl only when interaction is present. 
Assuming Co < z o, such smaller initial population imbal- 
ance leads to the same effect as a larger (smaller) tunnel- 
ing Cl > Cl (Cl < Cl), for repulsive (attractive) interaction 
17 > (f7 < 0). 

The scaling property is illustrated in Fig. [2j In (a), 
for a constant interaction N^U/Cl = 4.1, the solution for 
z^(0)/N^ = 1 is self-trapped, whereas for z^(0)/N^ = 
1/2, it oscillates. The scaled parameters U , CI allow us to 
recover the self-trapped behavior at a scaled population 
imbalance. In particular, the scaling properties allow us 
to understand why a small repulsive interaction leads to 
slower oscillations for 2^(0) = 1, while it leads to faster 
oscillations when z(0) < 1, as shown in Fig. [2^b). 




FIG. 2: (color online) Population imbalance scaling. (a) 
For NU/Cl = 4.1, the initial preparation z^(0)/N^ = 1 leads 
to self-trapping (black solid line), while for z^(0)/N^ = 0.5, 
oscillations emerge (blue dotted line). According to the 
scaling property Eq. (14), the parameters N^U /CL = 8.2, 
Cl = 0.13813, z^(0)/N^ = 0.5 recover the qualitative behav- 
ior of the solution z^(0)/N^ = 1 (red dashed line), (b) For 
z^(0)/Nijj = 1, increasing the interaction from N^U/Cl = 
(solid blue) to N^U/Cl = 1 (dotted blue) leads to a decrease 
of the observed frequency, whereas for z^(0)/N^ = 0.5, the 
opposite behavior is observed, which can be understood via 
the scaling property. 



We will assume r = 1, i.e. perfect isospecificity, unless in- 
dicated explicitly. A realistic physical implementation of 
the model consists in the second species being realized by 
a different hyperfine- state, i.e. one can convert particles 
between the species at will, by applying suitable Rabi 
pulses on one or both wells. We are interested in the 
dynamics of the spatial population imbalance when the 
second species (j) is populated at the expense of ip, i.e. for 
a constant total particle number N = + and con- 
stant initial spatial population imbalance z tot = z^ + z<j>. 



A. Equations of motion and conserved quantities 



III. BI-SPECIES DIMER 

Having established the behavior of a single-component 
BEC in the double- well, we now add a second component 
with very similar physical properties. The dynamics of 
the second species, which we will refer to as (/>, is governed 
by the same physical parameters as the first species 
the tunneling rates are equal, CL := Cl^ = f^, and the 
intra- and inter-species interaction strengths fulfill 

rU^rf = rU^ rU. (15) 



In analogy to Eq. ([T]), the coupled discrete Gross- 
Pitaevskii-equations for the bi-species system read 

i^j = ^(l^| 2 +r|0i| 2 )^--"^3-i, (16a) 
i^k = U(r\^ k \ 2 + \(f) k \ 2 )0 k - ^03- fc , (16b) 

where k,j = 1,2 and we assume r = 1. The coupling 
between the two species occurs via the non-linear inter- 
action term containing U. In the following, the total 



5 



particle number of each species is conserved, 
IV-iWI 2 + \Mt)\ 2 = Nip, |0i(t)| 2 + \4> 2 {t)\ 2 = Nt, (17) 

for all times. The two Bloch-vectors and that de- 
scribe the two species, defined analogously to Eq. ([2]), 
evolve on spheres of radius and Nj,, respectively. 



The coupled equations of motion Eq. (16) can be re- 



formulated for the Bloch-vectors and to give 



dt 



-U yp (z^ + 2 
Q Zf3 +U xp (z^ 



(18) 



where /3 = -0, 0, while the total energy 



iftot = 



ft 



r^tot- 



(19) 



is conserved. Since the total Hamiltonian of the system 
depends only on the spatial population imbalance, z to tj 
and on the sum of the ^-components of the Bloch-vectors, 
Xtot? ^ assumes the same form as for a single species, see 
Eq. (|6|. Consequently, the equations of motion for the 
total Bloch- vector v tot , 



Vtot = V^j + Vf, 



(20) 



possess the same form as for the Bloch-vector of a single 
species, Eq. ([5| [15 . The length of v tot is a constant 
of motion, which is a consequence of the isospecificity of 
the dynamics, i.e. of U^^ = = and = Q^. 
Note that while the dynamics of the total population 
imbalance is integrable, this is not necessarily the case for 
the dynamics of the individual species, i.e. the difference 
between the imbalances of the two species, v re \ = — 
#0, eludes a closed solution [21]. When the interactions 
are not isospecific (r ^ 1, U^^ ^ U<j>^), |#tot| is n °t 
conserved anymore. Similarly to the triple- well system 
[36H38] , chaos can then emerge in the two-species double- 
well- system [El EE] . 

The dynamics of the mixture can be computed in two 
different ways. We can either integrate the equations of 
motion of the individual species, Eq. (18), or we can ex- 



press the dynamics in terms of a single effective species, 
as suggested by the total Hamiltonian Eq. (19). In or- 
der to illustrate the underlying physical mechanism, we 
will start with the former and return to the latter in the 
following sections. 



B. Dynamical effects of a second species 



The integration of the equations of motion, Eq. (18), 



allows us to directly compare situations with one and 
two species. The dynamics of the spatial population im- 
balance z tot (t) of a single-species and a fully isospecific 
bi-species BEC is shown in Fig. [3] (a,b), the behavior of 



weakly an-isospecific mixtures can be observed in (c,d), 
we will discuss the relaxation from isospecificity in Sec- 



tion HIE below. The initial spatial population imbalance 
is set to Ztot(0)/N = 0.8, while we vary the population 
of the second species and the interaction U (see also the 
sketches above the graphs). We keep the total population 
N constant and an initially vanishing phase between the 
wells, for both species, see also Fig. [4j We assume that 
the two species correspond to different ground states of 
the same atoms and that a Rabi-pulse is applied on the 
right well, which transfers population from the first to 
the second species. In Fig.[3ja), no Rabi pulse is applied 
and only one species is present in the system, N^/N = 0. 
In (b) , a 7r pulse transfers all the right- well population to 
the second species, N^/N = 0.1. For intermediate val- 
ues of NU/tt (1 < \NU/tt\ < 6), populating the second 
species changes the dynamics of the spatial population 
imbalance z tot (t) strongly, the location of the oscillation 
regime on the NU /Q-axis is shifted to smaller values 
for increasing population of the second species. Given 
a fixed interaction parameter NU/Q, the behavior can 
therefore qualitatively differ for the different populations: 
In panel (el), for NU/Q = 4, the single-species BEC is 
clearly in the oscillation regime, while for N^/N = 0.1, 
no Josephson-like oscillations are observed anymore, the 
system is self-trapped. For NU/Vt = —2, the opposite 
behavior is observed: While for N^/N = 0.1, the system 
oscillates, it is self-trapped for the other scenario. 

We can understand the impact of the second species 
from the structure of the Hamiltonian Eq. (19): Let us 
first consider repulsive interactions, i.e. U > 0. Since 
we assumed x to t(0) > 0, the phase-coherence part of the 
Hamiltonian, — Qx tot (0)/2, is initially negative, which fa- 
vors oscillations. The total energy, H tot = L r z to t(0) 2 /4 — 
^#tot(0)/2, is reduced - for | if tot I < ^/2, the system 
is in the oscillation regime, since z tot = can then be 
attained energetically. In other words, initial coherence 
between the two wells with x to t(0) > favors oscilla- 
tions. Populating the second species in one well partially 
destroys the coherence, which increases the total energy 
and pushes the system towards self-trapping. 

Attractive interaction (U < 0) leads to the opposite 
effect: The interaction term in the Hamiltonian is nega- 
tive, and the phase-coherence part — Qx tot (0)/2 is benefi- 
cial for self-trapping. Destroying the phase coherence by 
populating a second species is beneficial for Josephson 
oscillations. In general, the competition of interaction 
and phase-coherence in the Hamiltonian for the initial 
states in (a) breaks the symmetry U — U present in 
panel (b). 

Summarizing, the behavior of the system can be 
changed from self-trapping to Josephson oscillations by 
the population of a second species, without changing 
£tot(0)/iV~, U or Q. As can be seen, the behavior of the 
two-species mixture and of the single-species condensate 
are similar - although subject to a shift on the NU/Vt- 
axis. 
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FIG. 3: (color online) Total population imbalance dynamics for (a) one species, (b) two perfectly isospecific species, (c,d) 
two slightly an-isospecific species with a ratio of the inter-to-intra-species interaction strength r = 1.1 (c), and r = 0.95 (d), 
i.e. = rll, see Eq. |l5| . The initial imbalance is z tot (0)/N — 0.8 always, and only the population of the right well is 

modified. In (a), the second species is absent, in (b,c,d), the population in the right well fully appertains to the second species 
(see sketches above (a,b,c,d)). The panels (el)-(e4) on the right show examples for the four situations and selected values of 
NU/Q (dot-dashed red: (a), dashed blue: (b), dotted black: (c), solid green: (d)). 



C. Effective tunneling coupling 

The effect of a second species on the mean field dynam- 
ics can be modeled by a single species system with an 
effective tunneling coupling Q. To derive this in analogy 
to Section II B we consider the dynamics of the spatial 
population imbalance z tot . Similarly to Eq. (J8|, we find a 
second-order differential equation with parameters that 
depend on the initial conditions: 



Ztot 



(21) 



Ztot "TT^ot - ^tot(O) 2 ) + fi (fi + Ux tat (0)) 



where x to t(0) and z to t(0) are defined through Eq. (120 
aty(0) can be related to z^(0) and i^(0) through Eq.T9 
which is valid analogously for £^(0), 2</>(0) and i^(D). 
We assume in the following i^(0) = i</>(0) = and 
x^(0), x ^(0) > 0. By comparison to Eq. (J8|, we find 
that the spatial population imbalance for the two-species 
condensate, z to t(t), exhibits the same time-dependence 
as the imbalance of a single-species BEC that is pre- 
pared with the same initial imbalance, z^(0) = ztot(0), 
but which experiences an effective tunneling coupling 
where 



O = 



fi (fi + U (y/N* - 2^(0)2 + - ^(O)^)) 



u 2 (N 2 -zt t(oy 



1/2 



u 

2~ 



N*-ztot(vy 



(22) 



is found by solving the resulting quadratic equation. For- where z^(t) is the solution to Eq. ffil. In other words, 

mally, z tot (t) satisfies the effect of breaking the phase coherence due to the 

z tot (jv^,jv , z^(o), 2 (o), £/,n,t) 

= z^(N^ + N<f,, z^(0) + Ztf,(0),U, O, t), (23) 
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addition of a second species is equivalent to a change 
of the tunneling coupling to Cl. Consistently with the 
findings discussed by means of Fig. [3| Eq. ([22} shows that 
Q < Q for repulsive interactions U > 0, while Q > Q for 
attractive U < 0. 



D. Effective particle number 

The incorporation of the effect of the second species in 
the modified tunneling coupling parameter Q in Eq. ( 22 ) 



relies on the exact form of the two- mode differential equa- 
tions Eqs. (|8| and (21). A more general way to under- 
stand the dynamics of z t ot? which will also be applicable 
in larger systems, can be found by considering the emerg- 
ing total Bloch-vector. In the bi-species case, the length 
of the Bloch-vector of either species reflects the respec- 
tive population, \v^\ = iV^, |^| = N^. The total number 
of particles in the bi-species case is the sum of the lengths 
of the two Bloch-vectors, \v^\ + = + = TV, 
while the length of the total Bloch-vector v tot , defined in 
Eq. (20), is, in general, smaller than N. 



The dynamics of a bi-species condensate of a total par- 
ticle number N is thus described by the Bloch-vector iT t ot 
that evolves on a sphere of radius 7V e ff, 



N e f{ = \v to t\ 



(24) 



i.e. the system behaves like a single-species system with 
a reduced total population iV e fj = |i?tot|- 

Formally, the spatial population imbalance in the bi- 
species case fulfills, 



z tot 

= z^Nes^^Q) + Z<f> (0),n,U,t), 



(25) 



where iV" e ff is given in Eq. (24), which should be com- 



pared with Eq. (23), where we found an effective tunnel- 



ing coupling Q. Applying the scaling of the Bloch- sphere 
radius, Eq. (11), and of the initial population imbalance, 



Eq. (13), we can reconcile the two perspectives (Eqs. (23) 
and (|25|)) with each other, 

H (^eff,^(0) + ^(0),O,C/,t) 



^ ^ (n, (^(o) + z (o)) , n, u, tj , 



(26) 



where ft is given by Eq. (14b). Inserting the parameters 



confirms that the effective particle number iV e fj indeed 
leads to an effective frequency Q with Q = consistent 
with Eq. (l22|). 



The geometric coupling of the Bloch-vectors allows 
us to understand the effect of concerting or obstruct- 
ing relative phases. When the relative phases of both 
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FIG. 4: (color online) Coupling of the Bloch-vectors that 
describe the two species. Left panel: (blue) and 
(red) couple to a total Bloch-vector vtot (black, dotted), as 
compared to a single-species Bloch-vector v\ (black, solid) 
of length \v^\ + \v$\ — + N<f, = N, with identical z- 
component zo. The length of vtot (dashed green Bloch-sphere) 
is reduced in comparison to the single-species case (thin solid 
black Bloch-sphere). Right panel: Depending on the relative 
orientation of and v^, the resulting effective vector length 
I #tot I varies. The configurations (a,b) correspond to the re- 
spective physical situations shown in Fig. [3] Only one species 
is present in (a); for (b) \v<f,\/N = 0.1 and and are anti- 
parallel, while in the intermediate case (c), \v^\/N — 0.05 
(not shown in Fig. 



species and the relative population imbalances are equal, 

0^(0) = zi>(P)/ N 1> = ^(O)/JV0, the dynamics of 

the two species are concerted, and the collective motion 
corresponds to the one of a single species with the same 
phase difference. Geometrically speaking, the Bloch- 
vectors and are parallel, and the mixture behaves 
like a single species with N = + N^. Likewise, 
z^/Nfj) = z^/N^ implies Cl = Q in Eq. (22). 

When 0^(0) ± 0^,(0) or z^(0)/N^ + ^(O)/7V the 



different initial phases or imbalances work against each 
other, and the motion of ip is then not concerted with 
the motion of </>. The dynamics can then be described by 
a single species with smaller population. Geometrically 
speaking, depending on their relative orientation, the 
length of the total Bloch-vector can be close to + iV^, 
or be significantly smaller, as illustrated in Fig. [4j As an 
example, we return to the configurations considered in 
Fig. [3j for which the coupling of the vectors is depicted 
in the right panel of Fig. [4j In both (b) and (c) the 
second species is initially localized in the right well and, 
consequently, z^ = — TV^, the Bloch-vector points down. 
As the population of the second species is increased, 
is tilted towards the z-axis and the length of the total 
Bloch-vector decreases. 

In general, the dynamics of the system is invariant un- 
der the simultaneous application of a Rabi pulse on both 
wells, which geometrically corresponds to a splitting of 
the constant total Bloch-vector + = v tot into two 
new parts v f , and v*., with v, 



V, = Vtot. 



In particular, one does not alter the system behavior 
when applying a Rabi pulse if only one well is occupied 
initially such that ztot(0) = d=iV. No phase coherence 
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is broken by the Rabi pulse, since no phase relationship 
between condensate wavefunctions in the left and right 
well is present in the first place. 



E. Robustness against deviations from 
isospecificity 

Since conventional optical traps operate with linearly 
polarised light that is far detuned from atomic resonance 
frequencies, atoms in different hyperfine states typically 
experience the very same optical potential, which moti- 
vates the species-independent tunneling rate Q. However, 
the species may be subject to slightly different inter- and 
intra-species scattering lengths [19] [39j [40] , i.e. r may dif- 
fer from unity (see Eq. ([15]) ) . The question arises how 
deviations from the ideal parameters impact on the sys- 
tem dynamics, since the treatment above relies heavily 
on the conservation of |i?tot|- Relaxing the assumption 
r = 1, the Hamiltonian becomes explicitly dependent on 

Vrel = - V(f), 



U { (l + r) 2 (l-r). 
H = j (—2— z tot + — 



xtou (27) 



and the dynamics of the total Bloch vector also depends 
on the relative vector v re \, 



dv t < 



dt 



u 



§ ((r - 1) y re \z re i - (r + 1) y to tZ to t) 



2 ((r + 1) Xtot^tot - (r - 1) Xrel^rel) + ^tot | (28) 



The length of the total Bloch- vector, \v tot | is therefore not 
conserved anymore and the time-derivative of its square 
reads 



d I- |2 



2(1 - r)U(x^ - x^y^z^. 



(29) 



Thus, due to the different intra- and inter- species inter- 
action strengths, the relative Bloch vector v re \ perturbs 
the motion of the total Bloch vector £f tot , which experi- 
ences an effective particle number that varies in time, or, 
equivalently, a varying interaction strength. 

The influence of the relative Bloch vector on the dy- 
namics is naturally constrained by 



\N^-N^\ < \v tot \ <N 4 , + N i , = N, 



(30) 



i.e. in general, the greater the population difference be- 
tween the two species, the smaller is the influence of the 
deviation from iso-specificity. 

On the other hand, although the amplitude of the 
perturbation induced by the relative vector in Eq. (28) 
is constrained by |(r — l)E/iV/4|, a system prepared on 
the separatrix, i.e. at the fragile borderline between self- 
trapping and oscillating regime (see Fig. [I]), can still be 
affected qualitatively by very small deviations around 
r = 1, which is apparent for NU/Vt = 2.5 in Fig. |3je2): 
While the singular non-oscillating case is attained for 



r = 1, the system is self-trapped for r = 0.95 and oscil- 
lates for r = 1.1. On the other hand, when the trajectory 
is stable, and well away from the separatrix, i.e. deeply in 
the oscillating or in the self-trapping regime, the change 
in the dynamics due to the an-isospecificity is very small, 
as can be observed in Fig. [3jc) and (d) and in the panels 
(el,e3,e4), for [7 = 4,-2,-4. 

Slightly differing inter- and intra-species interaction 
strengths can lead to chaos at the borderline between 
self-trapping and oscillation regime [21]. The effect of 
exchange symmetry breaking, however, dominates in the 
regular regime and outweighs the influence of weak an- 
isospecificity. 



F. Experimental implementation 

An experiment that implements a controllable bi- 
species system can be realized with current technology. A 
binary system that naturally offers the parameters that 
we require is given by the := \F = l^mp = 1) and 
|0) := \F — l,rap = —1) components of 87 Rb. Since 
spin-changing collisions are weak [41 ], we can assume 
that the population of either species is conserved and 
that they can be simultaneously trapped in the same po- 
tential [19]. Moreover, the intra-species interactions are 
identical, and the ratio of inter- to intra-species interac- 
tion fulfills U^/U^ w 1.0093 [HE]. De P endin S on 
the potential difference between the wells at the moment 
of the cooling, any desired spatial population imbalance 
can be achieved [13]. An imbalance of the population of 
the two species can then be implemented selectively on 
one well by using the techniques utillized in [42]. Here, a 
strongly focussed laser beam shifts the hyperfine transi- 
tion in one well compared to the other. The site-selective 
transfer is then realized using a microwave pulse at the 
shifted frequency. 

Using well-established experimental techniques, the to- 
tal population imbalance can be read off in a destruc- 
tive way [13] ; and when combined with a spin-dependent 
push-out [42 , the dynamics of the individual species can 
be observed. 

In conclusion, the population of a second species 
changes the balance between phase-coherence and inter- 
action energy. The new initial state of the system may 
lie in a possibly different phase-space. For very small 
(NU <C ft) and very large (NU ^> Q) interactions, the 
qualitative behavior is not changed dramatically, and ex- 
change effects are small: For very weak interactions (in 
the Rabi regime), each particle tunnels individually, and 
no collective effects emerge, whereas for very strong in- 
teractions, the inter-particle interaction outweights the 
exchange effects. In typical experimental regimes that 
lie between these extremes [13], however, the breaking 
of the exchange symmetry can have a strong qualitative 
impact on the system dynamics. 
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IV. MANY-BODY TREATMENT: 
SPIN-COHERENT STATES AND BEYOND 

The difference between single- and bi-species conden- 
sates can be understood by the competition of phase- 
coherence and interaction energy. To assess whether the 
physical arguments can be extended beyond situations 
that are well described by mean-field theory, we proceed 
to an analysis using second quantization. 



A. Two-well Bose-Hubbard model 



The Gross-Pitaevskii equation Eq. (16) is the classi- 



cal limit of the Heisenberg equations of motion retrieved 
from the Bose-Hubbard Hamiltonian, 



where a = <p. The operators Jj^ a fulfill the usual SU(2) 
commutation relations, and the Hamiltonian Eq. ( [31] ) can 
be re-written as 



U 



(33) 



where a constant summand has been omitted here with 



respect to Eq. (31). From Eq. (33), it is immediate that 
a coupled basis, i.e. the eigenstates of = Jk^ + Jk,<f> 
allows us to re-write the Hamiltonian in the Lipkin- 
Meshkov-Glick form [43H47] 



^BH = UJ Z 



(34) 



n 

'2 



( a i,^2,^ + a I,0 a 2,0 + H ' C -) > ( 31 ) 



where hk, a denotes the number operator, and at (hj a ) 
creates (annihilates) a particle of species a(= ip, (j>) in 
site j. Further below, we shall present the results of 
a numerical simulation of the dynamics induced by this 
two-species Hamiltonian. In order to investigate the re- 
duction to single-species dynamics, it will be useful to in- 
troduce an angular momentum representation, which can 
be done elegantly for each species using the Schwinger- 
boson model: 



equivalent to the one for a single species in the Bose- 
Hubbard dimer, in direct analogy to Eq. (19) - note that 
the z-projection m is related to z tot via 2m = z tot: which 



leads to different pre-factors in the consistent Eqs. (34) 
and (19). At this stage, the problem seems to reduce 



to finding the appropriate single-species configuration, 
just like in the mean-field representation. In the discrete 
case, however, the quality of the single-species reduction 
depends heavily on the initially chosen configuration, as 
we shall see. 



B. Bi-species spin-coherent states 



(^2,a^l,a ^l,c*^2,c*) 



(32a) j n foe mean-field approximation, binomial counting 

statistics of the many-particle state is assumed [48 , such 

(32b) that the treatment of the last Section is expected to be a 
good approximation to the spatial population imbalance 
exhibited by initial states that are spin-coherent, i.e. of 

(32c) the form 




where the number of particles TV is large and |vac) de- 
notes the zero- atom state. The expectation values of this 
state match precisely the population imbalances and 
z^ and the relative phases 0^ and 0$. Note that the rel- 
ative phase between the species is unobservable in our 
setting. 

We compare the time-evolution of the population im- 
balance of a state of the above form Eq. (35) with the 



classical description based on Eq. (18), in Fig. [5] The 
Hamiltonian Eq. ( 34 ) is integrated in a numerically exact 



way, without performing any further approximation. The 
initial spatial imbalance is z tot (0)/N = 0.9, the number 
of particles is N = 180, and the interaction parameter 
fulfills U/Cl = 1/60. We compare (a) the single-species 
situation to (b) the bi-species scenario. For the latter, 
we assume that a Rabi 7r-pulse on the right well leaves 
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FIG. 5: (color online) Population imbalance: Mean-field 
treatment and many-body dynamics. The initial population 
imbalance is z to t(0) = 0.9iV, U/Q = 1/60, N = 180. We com- 
pare a spin-coherent many-particle state of the form Eq. ( [35] ) , 
for (a) = 0, = N (all particles are of the same species) 
and (b) N^/N = 0.05, N^/N = 0.95 (particles in the differ- 
ent wells are of different species) to the classical calculation 
based on Eq. ( |18[ ). Black solid lines denote the classical cal- 
culation, red dashed lines the exact quantum bi-species cal- 
culation. The single-species description based on the effective 
tunneling constant Q is shown in blue solid, the one based 
on the effective total particle number iV e ff = 162 is shown as 
black dotted line. 



all particles in the second internal state 0, such that 
N^/N = 0.95 and N^/N = 0.05 (see inserted sketches). 
The classical calculation (black solid line) exhibits oscil- 
lations in (a) and self-trapping in (b) (note the changed 
scale for the ordinate). The quantum many-body cal- 
culation features the same qualitative behavior, but the 
oscillations decay due to quantum fluctuations [35] . 

The single-species descriptions based on the effective 

(blue solid line), or employ- 



coupling Q given in Eq. (22 



ing the effective particle number N e ff given in Eq. (24) 



(black dotted line), reproduce the dynamics well, and also 
feature the decay of oscillations. In particular, switching 
the population of the right well to the second species by a 
7r Rabi pulse clearly switches the behavior of the system 
from Josephson-oscillations to self-trapping - just like in 
Fig. [3] The decay of the oscillations for the description 
based on the effective particle number N e ^ = 162 (black 
dotted) is slower than for the two-species scenario (red 
dashed) and for the description via the effective tunnel- 
ing coupling Q (blue solid), since all N e Q particles start 
in the left well, which reduces quantum fluctuations. 



C. Bi-species Fock-states 

Since relative phase and particle number difference 
are conjugated variables, the phase relation between two 
wells populated by the same species can also be broken 
by preparing a Fock-state with a well-specified number of 
particles in each well. In the semiclassical picture, a Fock- 
state corresponds to a distribution on the Bloch-sphere 
with well-specified z-component, but uncertain x- and 
^/-component, i.e. to a ring on the Bloch-sphere, or a ver- 
tical line in Fig. [I] The resulting distributions for single- 
species and bi-species Fock-states differ, and so does the 
dynamics, as we will see below. 



1. Single- species description 

We now derive an effective single-species description of 

and 



a bi-species Fock-state, analogous to Section [III D 
discuss the requirements for its application. The Hamil- 
tonian Eq. (34) suggests to work in the basis of coupled 



angular momenta, i.e. in eigenstates of J 2 , J 2 , J|, J z , 



\j,3ij>J<t>,™>) , 



(36) 



where jy, = and = are related to the 

number of particles of either species, and m is the total 
spatial particle number imbalance. The quantum number 
j as well as the total particle number of either species, 
N^jj and TV^, are constants of motion - these conserved 
quantities are the analogous quantities to |i?tot|> Wip\, an d 
\V(p\ in the classical description. 

An initial Fock-state with a well-defined num- 

ber of particles of either species in each well is, in gen- 
eral, not an Eigenstate of J 2 of the form Eq. (36), but 



an eigenstate of the uncoupled angular momentum op- 
erators J 2 , J z ,ip-, «/|, Jz,(f>- It can, however, be written as 
a superposition of total angular momentum eigenstates 
given in Eq. (36), 



I^f) = b'^,^,j ,m ) = c 3 \3>U>M> m ) >( 37 ) 

3 = \U~U\ 

where m = + and Cj is the Clebsch-Gordan coef- 
ficient 



Cj = (jJi>Jcp,m\j<ip,rn^J ( f ) ,m ( f ) ). 



(38) 



The single-species description with a fixed effective to- 
tal particle number corresponds to the approximation 

I^f) = |j^,m^, j ,m ) « \ jeffJ4>U> rn 4> + m 1>) > ( 39 ) 

where j e ^ is the expectation value of the total angular 
momentum, 

ieffOeff + l) = <* F | J 2 |*F> (40) 

= U(U + !) + U(U + !) + 2^-^,7710, 
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FIG. 6: (color online) Time-evolution of Fock-states for N = 120, U/Q = 0.021. (a,b) Counting statistics for single- and 
bi-species Fock-state with initially 40 (80) particles in the left (right) well. Horizontal axis: Time t. Vertical axis: Number 
of particles in the left well. The color-code indicates the probability, (a) Bi-species Fock-state Eq. (43) with = = 
60, Ncf> = 2j(j> — 60, — 10, m<f, — —j^ — —30, m to t — —20. (b) Single-species state Eq. (44) with — 120, — —20. (c) 
Single-species approximation Eq. ( [45] ) for the bi-species state (a): N e s = 70, m = —20. (d) Average population imbalance, for 
the bi-species Fock-state (black dotted), single-species Fock-state (blue solid), the single-species approximation (red dashed) 
and the weakly an-isospecific case r = 1.1 (brown dot-dashed). 



which allows an interpretation as effective particle num- 
ber, iVeff/2 = jeff- The classical limit of Eq. (40) is in- 



deed Eq. (24), which neglects linear terms with respect 
to Eq. Q40[). 

The approximation Eq. (39) is reasonable only for 



states |^f) that possess a well-defined total angular mo- 
mentum, i.e. when 



A(J 2 

(P) 



«1, 



(41) 



such that the distribution of Clebsch-Gordan coefficients 



cj that appears in Eq. (37) is narrow and peaked at j e ff- 



The variance of the expectation value of J 2 amounts to 



" 3% 



)(jcj> + j% - ml) - Sm^m^, 



i.e. the uncertainty in the total angular momentum de- 
pends on the projection ra^, of the two angular mo- 
menta. Therefore, depending on the physical situation, 
the relation Eq. (39) is exact, approximate, or unsuitable 



in the limit of many particles: 



• Only one well occupied. We then have 
jtjj , ^0 = j(f> and find the equality 



\U,m^ fyjfamt = 3ct>) 



(42) 



Since all particles are prepared in the same well, 
the above state is fully exchange-symmetric and be- 
haves like a Fock-state of a single species. 

• One species localized in one well. If all par- 
ticles of one species are localized in one well while 
the particles of the other species remain distributed 
among the two wells, i.e. without restrictions of 
generality = j^, the relative variance Eq. (41) 
becomes small in the limit of large total particle 
numbers - we assume that N is increased for a con- 
stant ratio of species populations and constant pop- 
ulation imbalances (constant ra^/j^, m^/j^, and 
j(f)/jip)' Then, Eq. ( [39]) is approximately valid. The 
scaling property Eq. (|23|) , which was shown to hold 
for spin-coherent states, is here, however, not ap- 



plicable anymore: In the derivation Eq. (22) of ^, 
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we assumed a fixed value of 6^ and 0^. Since the 
initial phase relationship is uncertain - a Fock-state 
is prepared - no effective tunneling Cl emerges. 

• General state. If both wells are occupied by par- 
ticles of either species, we have \m^\ < and 
\ m 4>\ < 3<t>- The distribution of Clebsch-Gordan co- 
efficients in Eq. (37) remains broad for increased to- 



tal particle number, and, in particular, the distribu- 
tion of the Cj is not necessarily peaked at j e ff- The 
relative variance Eq. (41) does not vanish for large 



particle numbers, and the approximation Eq. (39) 
is unsuitable. 

We exemplify the validity of the approximation 
Eq. (39) for a two-species Fock-state in Fig.|6| for U/Q = 
0.021, N — 120. Panels (a-c) show the probability for a 
certain number of particles in the left well, and panel (d) 
compares the average population imbalance obtained by 
the single- and two-species numerical solution and by the 
approximation. In (a), the initial state is 

|*f,2> = \U = 3<W = 10^V = 30, m = -30) , (43) 

which leads to a very distinct counting statistics when 
compared to a single-species state with the same initial 
spatial population imbalance in (b), 



|^F,i) = |^=60,m^ = -20). 



(44) 



In particular, the probability to find almost all or almost 
no particle in one mode is clearly enhanced for |^f,i) 
with respect to |^f,2), well in accordance with our in- 
tuition for bosons. The typical features of the count- 
ing statistics of |^f,2) are well reproduced by the single- 
species approximation based on Eq. (|39|), 



(45) 



|^F,l,appr.) b'eff = 35, 771 = -20) , 



where the vertical axis of panel (c) is displaced for direct 
comparison to (a) and (b). Not only is the suppression of 
events with many particles in the left or right well and the 
overall structures of the counting statistics reproduced, 
also the characteristic oscillation of the center-of-mass 
agrees very well, as can be observed in (d). The more 
fine-grained interference pattern in (c) with respect to (a) 
can be understood from the unique contribution of one 
term in the sum Eq. (37), i.e. for |^f,2), many-particle 
interference is averagedout. 



2. Semiclassical description 

An intuitive argument for the validity of the single- 
species approximation Eq. (39) for = ±7'^ can be 
obtained in a semi-classical picture. A Fock-state cor- 
responds to a Bloch- vector with fixed z-component, but 
uncertain x- and ^/-component - the phase between the 
wells is maximally uncertain, since the particle number 
in each well is fixed [48-51 . A single-species Fock-state 




°\/N °- 5 



(b) 1.0 

0.8 

z 0.6 

N 0.4 

0.2 

0.0 
-1.0 



Zo/N 


-£jr^^ 


**» 

1 \ 


** 

*** 1 



-0.5 



0.0 /A7 0.5 

x/N 



1.0 



FIG. 7: (color online) Coupling of Bloch- vectors with speci- 
fied z-component, but random phase. We show a cut through 
the Bloch-sphere, for y — 0. Since the relative phases be- 
tween the wells are unknown for either species, the length of 
the emerging vector vtot is also uncertain. In (a), and 
couple to a total vector v to t with a length ranging between the 
small and the large dashed-green Bloch-sphere, since only the 
^-component, but not the x— and y— component of the two 
species are fixed. The vector is fixed in the representation 
for convenience. The effective total particle number is not 
well defined, but follows a distribution that corresponds to 
the classical limit of the Clebsch-Gordan coefficients. In (b), 
the x- and y- components of the 0-species vanishes, which also 
fixes the length of v to t- In the corresponding many-body sit- 
uation, the approximation Eq. (39) is valid for large particle 
numbers. 



is thus described by a vector of constant length that lies 
on a cone around the z-axis. A bi-species Fock-state cor- 
responds then to two coupled vectors, and of fixed 
length and fixed z-component, but uncertain x- and y- 
components. The z-component of the coupled, total vec- 
tor is well-specified, z tot = + z^. In general, however, 
the length of the total Bloch-vector is not fixed, since the 
relative orientation of and is uncertain. The bi- 
species system behaves like a single-species system with a 
total particle number that is not fixed, but which follows 
a certain probability distribution, P(\v tot \), which is pre- 
cisely the classical analogous [52] to the Clebsch-Gordan 
expansion Eq. (37). This situation is also depicted in 
Fig. 0(a). 



When one of the vectors, say points to the z- 
direction, it then exhibits vanishing uncertainty in the x- 
and ^-components - the cone on which it rotates shrinks 
to a point. Consequently, also v tot moves on a cone and 
has a constant, well-defined length, as can be observed 
inFig.^b). 
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3. Robustness against an-iso specificity 



A. Many-site Bose-Hubbard model 



Again, the question arises how robust the dynamics 
are with respect to weakly an-isospecific interactions, 
i.e. when r ^ 1 and thus U^^ ^ U^^ (see Eq. (15)). Due 



to the Heisenberg uncertainty relation, and due to z and 6 
being conjugate variables, no quantum state can be pre- 
pared precisely on the separatrix between self-trapping 
and oscillation regime. Thus, the pathological case shown 
in Fig. [3^e2) has no analogy in the quantum treatment. 

Since Fock-states can be described semiclassically as a 
distribution of classical solutions with fixed population 
imbalance z and unknown phase the overall behavior 
of the system is seldom clearly in the self-trapping or in 
the oscillation regime (consider a distribution shaped as 
a vertical line in Fig. IT]), and thus widely unaffected by 
small changes in r. We show the average spatial imbal- 
ance for r = 1.1 in Fig. [6] (d), it remains close to the 
behavior of the system with r = 1 (a). 

Summarizing, just like for the mean-field treatment, 
the impact of the exchange symmetry breaking can 
change the system behavior qualitatively, while inter- to 
intra-species interaction ratios that deviate little from 
unity do not jeopardize the general system behavior. 



V. EFFECTIVE SINGLE-SPECIES 
DESCRIPTION OF GENERAL MULTI- SPECIES 
BOSE-HUBBARD MODELS 



In the two previous sections, we have shown that the 
total particle density in a mixture can be be expressed 
as the density of an effective single species, either exactly 
- in the classical treatment - or approximately - in the 
quantum many-body treatment. On the one hand, this 
allows us to understand the effect of a broken exchange 
symmetry, on the other hand, especially in the many- 
body treatment, the relationship allows to save signifi- 
cant computational resources [53] . It is therefore appeal- 
ing to expand the discussion to more general scenarios. 

Whereas the non-linear scaling of the initial popula- 
tion imbalance, Eq. (13), is a peculiarity of the two- mode 
model, we can generalize the argument of the effective 



particle number to an arbitrary number of species and 
wells. 



For this purpose, we consider a general L-site Bose- 
Hubbard Hamiltonian that describes s distinct species, 



BH,gen 



j>k \a=l / 

5> ) 



i=i 



=i 



j=l \a=l / 



(46) 



where ej is the on-site energy for the j-th site, is 
the tunneling coupling between site j and fc, U is the 
interaction strength. The creation operatore a^ a creates 
a particle of species a on site j. The Hamiltonian Eq. (46 ) 
is isospecific, just like Eq. (31). 



In analogy to the Schwinger-Boson algebra, which al- 
lows to rewrite the two-site Hamiltonian in terms of the 



three generators of SU(2) (see Eqs. (32), (33)), we can 
re- write Eq. (46) in terms of the L 2 — 1 generators of 



SU(L), for each species a. This approach has been ap- 
plied to the triple-well case [54-56 and can be readily 
generalized to the present L-site system with s species. 
The L — 1 diagonal generators of SU(L) can be chosen 

as 



£(fc) 



Ek 
p=l n P,a 



k • hk+i ,, 



fc + 1 



(47) 



where k G {1, ... ,L — 1}, while the L(L — 1) hopping 
operators read 



y-(Z,m) _ 



(48) 



(49) 



where I > m. For L = 2, we simply retrieve Eq. (32). 
For a given species a, the L operators hj ja can thus 

(k) 

be expressed by the L — 1 operators Z\ ' and the total 
particle number operator N a = Y^i=i 



N 



L-l A(fe) 



(50) 



k=j 



where we set = for convenience. 

The SU(L) commutation relations of the x2 ,m ^ 
and Ya' m ^ are immediately inherited by the coupled op- 
erators 



ZW = £2< fc) , (51) 

s s 



a=l 
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since - just like in the two-site case - operators related 
to different species always commute. 

The Hamiltonian Eq. (46) can then be expressed with 
these coupled operators, 



the latter fulfills 



H 



BH,gen 



j>k 

L-l / k 



(52) 



Ne 



k=i 




where N is the total particle number, and e = 
(^^ =1 e)/L the average onsite-energy. The constant 
terms proportional to N and N 2 can be neglected for the 
dynamics. Just like for L — 2, the many-species Hamil- 
tonian assumes the same form as for a single species. 
Again, any multi-species state can be expanded in single- 
species states: A many-particle state of N a particles is 
defined by L quantum numbers, i.e. a general Fock-state 
with a well-defined number of particles in each mode is 
described by 

|^F,gen) = ®a=l 1^,^,1,^,2, • • ^^L-l) , (53) 

where m a ^ is the quantum number corresponding to the 
Z a operator and N a is the total particle number of the 
species a. In analogy to Eq. ( |37| ), |^F,gen) can be ex- 
panded in a sum of total particlermmber states, 

I* F,gen) = (54) 
s s 

Yl CA/ tot l^tot, • • • , N 3 , ^2 m a> i, . . . , ^2 m ajL -i> , 



where the CN tot are the generalized Clebsch-Gordan co- 
efficients for SU(L) [57] . 



B. Classical limit 



The Hamiltonian Eq. (52) allows us to retrieve the 
Heisenberg equations of motion for the s(L 2 — 1) time- 
dependent operators in the Heisenberg picture, exploit- 
ing the commutation relations of the generators of SU(L). 
In the mean-field limit, these equations become discrete 
Gross-Pitaevskii equations for vectors v a on s (L 2 — 1)- 



dimensional Bloch-spheres, analogous to Eq. (18): 



dt 



Va = M\J2 



V(3 V a , 



(55) 



,/3=l 



where M is a matrix of dimension (L 2 — 1) x (L 2 — 1). 
Since M depends only on the total Bloch-vector 



^tot — Y2 Va ' 



(56) 



dt Vtot 



M(v t ot)v t0 f 



(57) 



That is to say, the inclusion of a new species with the 
same physical properties as the existing one leads, again, 
to a total Bloch-vector i? tot , whose length can be inter- 
preted as an effective particle number. 



VI. CONCLUSIONS 

Ultracold atoms in optical lattices do not only consti- 
tute ideal model systems to study the effects of tunable 
inter-particle interaction, e.g. via Feshbach resonances 
[58] , but also to investigate the consequences of a broken 
exchange symmetry or distinguishability of particles via 
the population of different hyperfine states. Thereby, the 
interplay of interaction and exchange symmetry - the two 
pertinent ingredients of quantum many-body dynamics - 
may be explored. We investigated a two-well system, in 
which the population of a second internal state effectively 
breaks the phase relationship between the macroscopic 
wavefunction in the two wells, and thus influences the 
dynamics of the system. The bosonic Josephson junc- 
tion exhibits symmetries that allow us to understand the 
dynamics either by an effective tunneling coupling, or by 
an effective total particle number. 

For double-Fock-states, for which no phase relation 
between the wells exists, exchanging particles of one 
species by another breaks the exchange symmetry of the 
total many-particle state. The effective particle num- 
ber is then, in general, not constant, but follows a 
probability distribution, which affects the many-particle 
interference-capability of the system, and which can be 
observed in the emerging counting statistics. 

We showed that controlling a degree of freedom that is 
invisible for the Hamiltonian of a system can effectively 
steer and influence the dynamics considerably by break- 
ing or restoring phase coherence. Here, the dynamics in 
the double-well can be switched between the oscillatory 
and the self-trapping regime, which may be used as an 
effective interaction-controlling parameter, e.g. when Fes- 
hbach resonances are unaccessible. While we have con- 
fined ourselves to the comparison of situations with con- 
stant species population, driving a system via breaking 
and restoring its phase coherence may allow for further 
control [59] . 

By relaxing the assumption of perfect isospecificity, we 
also showed that the breaking of the exchange symme- 
try has a much greater impact on the system dynam- 
ics than weakly differing inter- and intra-species interac- 
tions. Only in the fragile case of orbits near the transi- 
tion between dynamical regimes does a small change in 
inter- to intra-particle interaction strength jeopardize the 
single-species model. With the inclusion of more species 
simulations become more and more demanding, and a 



15 



single-species approximation can potentially offer a sig- 
nificant computational advantage. It remains to be stud- 
ied to which extent quantum phases of multicomponent 
bose gases can be described by methods analogous to 
those exposed here. 
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